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Abstract. Over the last decade, the Concentration, Asymmetry and Smooth- 
ness (CAS), as well as the M20 and GINI parameters have become popular to 
automatically classify distant galaxies in images. Ellipticals, spirals and irregu- 
lar galaxies all appear to occupy different regions of this parameter space. At 
the same time, the Source Extractor (SE) program has become the mainstay 
to produce object catalogs from large image surveys. A logical next step would 
be to incorporate the structural parameters into the Source Extractor software. 
There are however several problems that arise: 1) the CAS parameters are fits 
to the images and Source Extractor eschews fits in the interest of speed, 2) the 
definition of the structural parameters changed over time. Now that there is 
a clear and agreed-upon definition of the structural parameters, I am incorpo- 
rating computed versions in the Source Extractor code {v2.5). The fitted CAS 
parameters are available for the GOODS-N/S fields and I compare the com- 
puted structural parameters to those found by the previous fits. My goal is to 
expand the source structure information in Source Extractor catalogs in order 
to improve automatic identification of sources, specifically of distant galaxies. 

The computed parameters perform reasonably close to the fitted versions 
but noise appears in faint objects due to a lack of information. For a subset 
of objects, the asymmetry signal is outside the SE boundaries and Smoothness 
still fails to compute for many objects. Type classification based on the SE 
parameters still lacks resolving power. 



1. Galaxy Structural Parameters 

There are six parameters that are popular for galaxy classification: Concentra- 
tion (C), Asymmetry (A), Smoothness/Clumpiness (S), Moment of the top 20 
% pixels (M20), the GINI-parameter (G) and the ellipticity of the object (E): 
Concentration: Conselice (2003) defines Concentration as the 5 x logiQ{r^Q / r2Q) 
with rgo and r2o the circular radii encompassing 80% and 20% of the flux are 
computed by Source Extractor (FLUX_RADIUS with PHOT_FLUXFRAC at 
0.8 and 0.2). Different flux percentages are sometimes used, e.g., the SDSS uses 
f^o/f^o- Since the computation of radii is done in SE, one would expect an 
straightforward implementation of the Concentration parameter in SE. 
Asymmetry: Asymmetry is the absolute difference of the original object (1) 
with the same object, rotated by 180° (R), and divided by the total flux of 
the object: A = abs{I — R)/I- Asymmetry is normally fit to objects with 
the X and y values of the center of rotation as variables. The information - 
position and value of all the pixels belonging to an object- is available in the 
SE data-structure. Thus, implementation of a calculated -not fitted- version of 
asymmetry is possible. 
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Smoothness: Smoothness is similar to Asymmetry but one subtracts a smoothed 
version of the object from itself. This parameters definition has changed the 
most over time with both shape and size of the smoothing kernel changing sig- 
nificantly. At present the smoothing is done with a boxcar smooth using 0.1 
Petrosian Radius as the smoothing kernel. 

M20: Lotz et al. (2004) introduced the second-order moment of the brightest 
20% of the galaxys flux (M20 = logio (SMi/Mtot) , with Sfi < 0.2 ftot) compared 
to the total second moment of the object (Mtot = Sffi[(xi — Xc)^ -|- (yi — yc)^])- 
SE already works with second order moments but the implementation would 
need an ordered list of pixel values (fi) and positions {xi and Ui). 
GINI: The GINI parameter is the second parameter used by Lotz et al. (2004) 
and Abraham et al. (2003). It is the area between the distribution of pixel 
values and a uniform distribution of pixel values. Lotz et al. (2004) present 
a computationally cheaper version which uses a ordered list of pixel values: 

EUipticity: Ellipticity is defined by the Source Extractor manual as: E = 1 — 
b/a, where a is the major and b the minor axis of the object. SE computes the 
values of a and b from the second-order moments of the pixels. 

The above parameters, except Smoothness, are used in the Zurich Estimator 
of Structural Type (ZEST, Scarlata et al. 2007). They determined the principal 
components of the above parameter-space from SDSS and COSMOS data and 
published the eigenvectors. So given the above parameters, one could easily 
determine the ZEST galaxy type-classification. 



2. Implementation and Test Data 

The above parameters are implemented in Source Extractor^ v2.5 code (SE 
Bertin and Arnouts 1996; Holwerda 2005) with two general classes: object and 
FIELD versions. The object version is computed from the pixels that SE has 
assigned to the object. The FIELD version is computed within the 1.5 Petrosian 
radius aperture. 

The Great Observatories Origins Deep Survey (GOODS) provides us with 
2 deep mosaics made with HST/ACS in several filters. CAS parameters were 
determined by Chris Conselice and collaborators (sec also Bundy et al. 2005). 
The SE parameters and catalogs used by the GOODS team are also public^. I 
re-ran SE with the same parameters as the GOODS team but now with added 
structural parameters on the central fields of GOODS-S (22, 23, 24, 32, 33 
and 34 in z filter). A cross-correllation of the catalogs gives the control for 
SE implementation. The SE run took substantially longer with the structural 
parameters than without. FIELD parameters substantially slow down SE, object 
ones inflict no performance penalty. 



^http : / /terapix . iap . f r/ 

^http : //archive . stsci . edu/pub/hlsp/goods/catalog_r 1/ 
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Figure 1. The relations between SE values and those determined by C. 
J. Conselice (CAS), left panels: eoncentration indices from SE C8020 = 
^80/^20 a-nd C9050 = rga/r^Q, and the concentration value provided by C. 
Conselice. The Cso20 still suffers from an offset, middle panels asymmetry, 
both the object and computed with 1.5 Petrosian radii (FIELD), compared 
to the CAS value. The FIELD version can be computed for just any entry 
in the catalogs, resulting in a lot of noisy entries with A^O.5. right panel: 
the smoothness parameter only resulted in overlap with the CAS catalog for 
the FIELD version. Smoothness fails repeatedly in either catalog, making it 
unreliable for now. 



3. Performance of Structural Parameters 



Figure 1 shows the relation between the values from the CAS fits and the val- 
ues determined by the SE on GOODS-S. The relation between concentration 
parameters indicate that the GOODS CAS catalog simply uses C = rsQ/r2o- 
The Asymmetry parameter becomes very much noisier in the FIELD version. 
Smoothness measures fail often for the original CAS fit (S=0) but also for the SE 
implementation. Only the FIELD version of Smoothness produces some overlap. 

Figure 2 shows the difference between the SE determined value and those 
from the CAS catalog for Concentration, Asymmetry and Smoothness. Noise 
increases at lower fluxes, as can be expected with the dearth of information for 
faint objects. For some large objects, most of the Asymmetry information is 
outside the object as defined by SE. 

ZEST Type classification based on structural types seem possible as well as 
an estimate of Sersic index (n) from the ratio of Petrosian and effective radius 
(rso, Graham et al. 2005). However, the SE classification does not differentiate 
types well enough yet. 
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Figure 2. The difference between the CAS and the SE values as a function 
of flux. As expected, at lower flux, the parameters become erratic. Con- 
centration suffers from an unexplained offset. Asymmetry is larger for CAS 
some cases because some information is outside the group of pixels assigned 
to the object by SE. The Smoothness_FIELD parameter works, but only for 
a fraction of the objects in the catalog. 



4. Conclusions 

Based on the first application of the modified SE v2.5 on the GOODS-S fields, 
I conclude the following: 

• Galaxy structural parameters can be implemented in Source Extractor, 
opening up possibilities for structural classification in large-scale surveys. 

• Concentration values suffer from an unexplained offset (Fig. 1). 

• Faint emission at the edge of objects may dominate Asymmetry (Fig. 2). 

• Smoothness fails often in both the original catalog and the SE approxima- 
tion(Fig. 1) but agrees when computed (Fig. 2). 

• M20 and GINI are clearly defined and appear to work. 

• ZEST galaxy type based on SE is not robust yet (all objects are t~2). 
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